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We calculate the strong forni factors and coupling constants of D* DsK and D*DK vertices using 

the QCD sum rules technique. In each case we have considered two different cases for the off-shell 

particle in the vertex: the ligthest meson and one of the heavy mesons. The method gives the same 

coupling constant for each vertex. When the results for different vertices are compared, they show 

. . . that the SU(4) symmetry is broken by around 40%. 
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; I I The knowledge of the form factors in hadronic vertices is of crucial importance to estimate hadronic amplitudes 

Qh' when hadronic degrees of freedon are used. When all the particles in a hadronic vertex are on-mass-shell, the effective 
■^ ' fields of the hadrons describe pointlike physics. However, when at least one of the particles in the vertex is off-shell, 
the finite size effects of the hadrons become important. 

In this work we study the D*DsK and D*DK vertices, which are fundamental to the evaluation of the dissociation 

cross section of J/^ by kaons when using effective Lagrangians. The supression of charmonium production is one of 

the most tradicional signatures of the quark-gluon plasma (QGP) formation in relativistic heavy ion collisions 1]. The 

J> [ dissociation of charmoniums in the QGP due to color screening would lead to a reduction of its production in such 

C t collisions. However, using the charmonium suppression as a signature of QGP formation requires the understanding of 

^^ ' J/^ production and absortion mechanisms in hadronic matter, because this suppression may be due to the interactions 

^~r , with the comovers during the collision |^ |^ . 

i^i. ■ One of the approaches used to study the interaction of charmonium with the hadronic medium, mainly in the low 

\^ i energy region {^/s < lOGeV), is based on effective SU(4) Lagrangians 0, S B S • This technique, however, 
f^ ' requires the detailed knowledge of the form factors in the hadronic vertices. The calculated cross section may change 
by a factor of two if a soft, instead of a hard, form factor is used in the vertices containing charmed mesons. 

This situation gave us the motivation to start a program to calculate charmed form factors and coupling constants, 
using the QCD sum rules approach lOJ. We have been continuously working on this problem and computing different 
M-i, vertices [III [12, 13, 14, 15, 16, 17J. An interesting subproduct of such calculations, 12, 13, 14, 17], was the undcr- 
Zh ■ standing of the behavior of the off-shell particle probing of the vertex: heavier particles resolves better the structure 
... of the vertex, while lighter particles are more suitable for measuring its size. This conclusion is also supported in the 
^ ■ present work. 

k>( I As a part of this project we evaluate, in the present calculation, the form factors in the vertices D*DsK and D*DK, 

5_j ■ and compare the results with the predictions from the exact SU(4) symmetr y [9 1 . 

d ' Following the QCDSR formahsm described in our previous works [lil[l2j[l3;[lJ>[la>[la[l3) we write the three-point 

correlation function associated with the D* D^K vertex, which is given by 

rf b,P') = j d^xd'v e^^'-^e-^(p'-P>y{Q\T{jf{x)j''\y)f'\Q)m (1) 

for K meson off-shell, where the interpolating currents are j^ — cry^d, j^ = is'j^d and j^" = icj^s, and 

ri'l^Hp,?') = I d'xd'y e^^-e-^^-^-y (0|r{jf (x)j^=^(y)j,^*^(0)}|0) (2) 






for Ds meson off-shell, with the interpolating currents j^ = vrj^^^s, j^' — icj^s, jj^ — u-ff^c, with u, d, s and c 
being the up, down, strange and charm quark fields respectively. In both cases, each one of these currents has the 
same quantum numbers as the corresponding mesons. 

Using the above currents to evaluate the correlation functions ^ and ||2Jl, the theoretical or QCD side is obtained. 
The framework to calculate the correlators in the QCD side is the Wilson operator product expansion (OPE). The 
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FIG. 1: Perturbative diagrams for K off-shell (left) and Ds off-shell (right) correponding to the D*DsK vertex. 



Cutkosky's rule allows us to obtain the double discontinuity of the correlation function for each one of the Dirac 
structures appearing in the correlation function. Calling pi the spectral density for the Dirac structure i, we can 
write the correlation function as a double dispersion relation over the virtualities p^ and p' , holding Q^ = —q^ fixed. 
Therefore, the amplitudes Ti are given by: 



1 /-So /-"o 

T,{p\p'\q^)^-\ ds 



du 



Pi{s,u,Q'^ 



is-p^){u-p'')' 



(3) 



where the spectral density pi{s, u, Q^) equals the double discontinuity of the amplitude ri{p'^,p' , Q^). The amplitudes 
receive contributions from all terms in the OPE. The leading contribution comes from the perturbative term, shown 
in Fig. n The phenomenological side of the sum rule, which is written in terms of the mesonic degrees of freedom, is 
parametrized in terms of the form factors, meson decay constants and meson masses. The QCD sum rule is obtained 
by matching both representations, using the universality principle. The matching is improved by performing a double 
Borel transform on both sides. 

The perturbative contribution for both Eqs. ^ and Q, written in terms of Eq.Q, is given by 



for K off-shell, and 



|pp [A {rn\ — mciris — 2k ■ p + p ■ p'^ + 2tt (^rn^ — fc • p')] 



27r\/A 

+ p'^ [B (m^ — rncirLs — 2k ■ p + p ■ p'^ +2it {—m^ + mciris + k ■ p)] } 



(4) 
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{g^ii' [tt {ms (s - ml) - iTfic [u - ml)) + 2D {mg ~ TOc)] 



+ {PiiP'u +p'f,Pi^) [Amc - Bmg + 2C{ms - m^)] 
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(5) 



for Ds off-shell. Here s = p , u = p , t 



A = A(s, i, u) == s^ + t^ + M^ _ 2st - 2su - 2tu, k ■ p 
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and A, B, C, D, E and F are functions of {s, i, u}, given by the following expressions: 
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TABLE I: Parameters used in the calculation of the QCD sum rule for the D* DgK vertex. All quantities are in GeV. 



for K off-shell, and 
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2^/i ' 



\k\ = /co, 
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for Ds off-shell. 

The phenomenological side of the vertex functions is obtained considering the contributions of the Ds and D* 
mesons to the matrix element in Eq. |^ and the D* and K mesons to the matrix element in Eq. ^. We introduce 
the meson decay constants /;<-, f^^ and fo*, which are defined by the following matrix elements: 



{^\r\K) = 
mu'\D*) 



mj^L 



kJk 



Wig + TUq ' 



'-D. 



D^ 



nic + TUs 



(6) 

(7) 
(8) 



where e^ is the polarization vector of the D* meson. 

In principle, we can work with any Dirac structure appearing in the amplitude in Eqs. 1^ and ^. However, there 
are some points that one must follow: (i) the chosen structure must also appear in the phenomenological side and 
(ii) the chosen structure must have a stability that guarantees a good match between the two sides of the sum rule. 
The structures that obey these two points are p'^, in the case K off-shell, and p'^p'y in the case Ds off-shell. The 
corresponding phenomenological amplitudes in these structures are 



rwp'^(p2^p'2, 
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fD'fDjKmo.mlml 
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(9) 



for the K off-shell, and 



T(^^^^'{p',p'\Q')^g\,%AQ 



{-2)ifD* fo, fnrriD* m% 



{trie + ms){p'^ - rn\i.){p' - mj^)iQ'^ + m% ) 



(10) 



for Ds off-shell. 

In the case of K off-shell the contribution of the quark condensate vanishes after the double Borel transform. In 
the case of the Ds off-shell, the quark condensate does not contribute to the chosen structure. 

To write the sum rules we equate each phenomenological amplitude in Eqs. (j ^ - (|l()|l . with the expression obtained 
by substituting the corresponding spectral density in Eqs. Q)-© into Eq. (jSJ. The matching between both sides is 



M^ and P'^ = -p 



a 



M'' 



improved by performing a double Borel transformation pL8l | in the variables P^ = —p^ 

We get then the final form of the sum rule, which allow us to obtain the form factors g£i*£, k{Q^) appearing in 
Eqs. iP)- ((Tn)l . where M stands for the off-shell meson. 



We use Borel masses satisfying the constraint M^ jM'^ = rl\l^|ra^^^, where mi„ and 



are the masses of the 



incoming and outcoming meson respectively. In the case of the K meson off-shell, this constraint gives AP /M'"^ = 



^IJr 



ijj. 



For the Ds meson off-shell, the relation should be M /M'' 



• /m\. However, the small value of the 



K mass spoils the stability of the Borel transformation. Thus, as is common in the literature, we change the K mass 
for the p mass. The resulting relation is then IvP /M'"^ — rrijj, /m? . 

The values of the parameters used in the calculation of the D*DsK vertex are despicted in TableQ] The continuum 
thresholds Sq and Uq, appearing in Eq. ((2Jl, are given by sq — {rrim + A^)^ and uq — {rriout + ^uY , where rriin and 
ruout are the masses of the incoming and outcoming mesons respectively. For the K off-shell we have min — rno^ and 
mout = rnu*, and for the Ds off-shell we have rriin = mu» and mout — rnx (see Fig. ^. 

Using As = A,j = 0.5 GeV for the continuum thresholds and fixing Q^ — 1 GeV , we found a good stability of the 
form factor Qj^Jjj ^, as a function of the Borel mass M^, in the interval 3 < M^ < 5 GeV^, as can be seen in Fig. |21 
In the case of the form factor g^,/^ ^ the interval for stability of the sum rule is 2 < M^ < 5 GeV^, as can be seen in 



Fig. 13 Fixing Ag = A„ = 0.5 GeV and M^ = 3 GeV^, we evaluate the momentum dependence of both form factors. 
The results are shown in Fig. ^ where the squares corresponds to the QdJ^ k^Q'^) form factor in the interval where 

the sum rule is valid. The triangles are the result of the sum rule for the g^j.'o k(Q^) form factor. In the case of the 
K meson ofF-shcll, our numerical results can be parametrized by an exponencial function (dotted line in Fig. Q: 



g'D,lKiQ')^^-»Se~^- 



(11) 



As in Ref . [l^l , we define the coupling constant as the value of the form factor at Q^ 
of the off-shell meson. For the K off-shell case the resulting coupling constant is: 



I'M' 



where ttim is the mass 



9d*d,k 



3.01. 



(12) 



In the case when the Ds meson is off-shell, our sum rule results can be parametrized by a monopole formula (solid 
line in Fig. 0J: 



9d*"d^k\Q 



9.01 



g2 + 6.86' 
giving the following coupling constant, obtained at the D^ pole: 



(13) 



0^-°=^ - 3 02 



(14) 



Comparing the results in Eas. H12l) and 114|l we see that the method used to extrapolate the QCDSR results in both 
cases, K and Ds off-shell, allows us to extract values for the coupling constant which are in very good agreement with 
each other. 

In order to study the dependence of these results with the continuum threshold, we vary As = A„ in the interval 
0.4 < As = A„ < 0.6 GeV, as can be seen in Fig|31 This procedure give us uncertainties in such a way that the final 

results for the coupfings in each case are: gj^Jjj ^ = 3.02 ± 0.15 and g]^/^ ^ = 3.03 ± 0.14. 

Now we study the D*DK vertex. The treatment is similar to the previus case. The correlation functions are 



rfHp^P') = fd'xd'y e^^'-e-»(P'-P>^(0|r{j^''(x)j^^y)j^\0)}|0) 



(15) 
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FIG. 2: Stability of goJo^KiQ = 1 GeV''), as a function of the Borel mass M^ 




FIG. 3: Stability of g^o-D.AQ = 1 GeV''), as a function of the Borel mass M^ 
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FIG. 4: gjj.jj jf (squares) and 315.0 ^ (triangles) form factors as a function of Q from the QCDSR calculation of this work. 
The solid (dotted) line corresponds to the monopole (exponential) parametrization of the QCDSR results for each case. 



for K meson off-shell, where the interpolating currents are j^ ° — cjf^s, j^ — iwy^s and j^ — icj^u, and 



for D meson off-shell, with the interpolating currents jf^ = ^T/jTss, jii " — cj^s, and j^ — iu'y^c. See Fig. El for 
understanding the pcrturbative contribution with these currents. 

For each correlation function, Eqs. H15() and H16|l . the corresponding perturbative spectral density which enters in 
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FIG. 5: Dependence of the form factor with the continuum threshold, for K and Ds off-shell cases. The dotted line correponds to 
As = A„ = 0.4 GeV, the solid line corresponds to As = A^ = 0.5 GeV and the dashed one corresponds to As = A^ = 0.6 GeV. 



Eq. © is: 
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for D off-shell, where k ■ p 
the D*D^K vertex, with 
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for D off-shell. 

The phenonienological side of the vertex functions are obtained by considering the contributions of D and 13* 
mesons to the matrix element in Eq. I|15() and D* and K mesons to the matrix element in Eq. 1)16(1 . We introduce the 
decay constants fo and /d* , which are defined by the following matrix elements: 
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FIG. 6: Perturbative diagrams for K off-shell (left) and D off-shell (right) correponding to the D*DK vertex. 
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TABLE II: Parameters used in the calculation of the QCD sum rule for the D*DK vertex. All quantities are in GeV. 



where e^ is the polarization vector of the D* meson. The fn decay constant was already defined in Eq. © . Again 
we have to choose one Dirac structure for each case in Eqs. (|17|I -H18 () . Following the points discussed before, the 
chosen Dirac structures are p' for the off-shell K, and p'^^p^ for the off-shell D. The corresponding phenomenological 
amplitudes in these structures are 



T^^^'\p',p'\Q')^9'^%AQ' 
for K off-shell, and 
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(22) 



for D off-shell. As in the case of the D*DsK vertex, the quark condensate does not contribute to the sum rule for 
these structures. 

The procedure to obtain the QCD sum rule is the same used in the case of the D*DsK vertex studied before. 
In this case we use the following relations between the Borel masses: M'^/M''^ = mjj/m'jj, for K off-shell and 
A'P /M'"^ — m?jjt /m? for D off-shell. The values of the parameters used in the calculation of the D*DK vertex are 
given in Table llll where we have used the relation fo* ~ fo'lD^/fo and the value of fos/fo from Ref. [22] in order 
to obtain the D* decay constant. 

Using As = A„ — 0.5 GeV for the continuum thresholds and fixing Q^ — 1 GeV^, we found a good stability of the 
sum rule for gjyJjjxy ^^ ^ function of the Borel mass M^, in the interval 2 < AP < 5 GeV^, as can be seen in Fig. |7| 
In the case of gjjJjjjf, the interval for stability is also 2 < M^ < 5 GeV^, as can be seen in Fig. |S| 

Fixing As — A^ — 0.5 GeV and M^ = 3 GeV^ in both cases, we calculate the momentum dependence of the form 
factors which are shown in Fig. O The squares corresponds to the gj^DfciQ^) form factor in the interval where the 

sum rule is valid. The triangles are the result of the sum rule for the g^tDxiQ^) form factor. In the case when the 
K meson is off-shell, our numerical results can be parametrized by an exponencial function (dashed curve in Fig. O: 



5D'^L(Q') = 2.69e- 
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The coupling constant was obtained as the value of the form factor at Q^ 
constant is 

9d'DK = 2.87. 



(23) 
-mj^. In this case the resulting coupling 
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In the case when the D meson is off-shell, the sum rule result is represented by the triangles in Fig. [^ and they can 
be parametrized by a monopolc formula (solid line in the figure): 
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FIG. 7: Stability of g^.j^j^iQ = 1 GeV"'), as a function of the Borel mass M^ 
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FIG. 8: Stability of QhioxiQ = 1 GeV^), as a function of the Borel mass A'P 



giving the following coupling constant, obtained at the D pole: 

0^^^ - 2 72 



(26) 



Studing the dependence of our results with the continuun threshold, for Ag^u varying in the interval 0.4 < As_„ < 
0.6 GeV, as can be seen in Fig. ^| we obtain the following values, with errors, for the couplings in each case: 
c,(^^^ = 2.87 ± 0.19 and g^o-DK = 2-72 ± 0.31. 
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FIG. 9: ffjj.jjx (squares) and gjjJj^j^ (triangles) form factors as a function of Q^ from the QCDSR calculation of this work. 
The dashed (solid) line corresponds to the exponential (monopole) parametrization of the QCDSR results for each case. 
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FIG. 10: Dependence of the form factor with the continuum threshold, for K and D off-shell cases. The dotted line correponds 
to As — A„ = 0.4 GeV, the solid corresponds to As = A„ — 0.5 GeV and the dashed corresponds to As = A„ = 0.6 GeV. 



Concluding, we have estudied the form factors and coupling constants of D*DsK and D*DK vertices in a QCD 
sum rule calculation. For each case we have considered two particles ofF-shell, the lightest and one of the heavy ones: 
the K and Dg mesons for the D*DsK vertex, and the K and D mesons for the D*DK vertex. In the two situations, 
the off-shell particles give compatible results for the coupling constant in each vertex. The results are: 



gO'DK 



3.02 ±0.14 
2.84 ±0.31 



(27) 
(28) 
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We can compare our result with the prediction of the exact SU(4) symmetry [aHISj which would give the following 
relation among these numbers Q: go'DsK = 9d*dk = 5. Eqs. (|77|l and 1^5)1 shows that the coupling constants in 
the vertices D*DsK and D*DK are consistent one with the other, but that they are relatively far from the value 
given by the SU(4) symmetry in the cited reference. Therefore, we conclude that the SU(4) symmetry is broken 
by approximately 40% in the calculation performed here. We can also extract the cutoff parameter, A, from the 
paramentrizations in Eqs. (fTT|l and f^ for K off-shell, Eq. (fU^ for Ds off-shell and Eq. ((JSJ for D off-shell. We get 
A « 2.07 GeV for the K meson off-sheU , A w 2.61 GeV for the D^ meson off-shell, and A w 2.51 GeV for the D 
meson off-shell. Comparing the values of the cutoffs, we see that the form factor is harder if the off-shell meson is 
heavier, implying that the size of the vertex depends on the mass of the exchanged meson: the heavier is the meson, 
the more as a point like particle is its behavior when probing the target, as observed in Refs. p^ lla. Il4. Il7l| . 
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